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ABSTRACT 

We have examined the effect of slow growth of a central black hole on spherical galaxies that obey 
Sersic or R 1 / m surface-brightness profiles. During such growth the actions of each stellar orbit are 
conserved, which allows us to compute the final distribution function if we assume that the initial 
distribution function is isotropic. We find that black-hole growth leads to a central cusp or “excess 
light”, in which the surface brightness varies with radius as R~ 13 (with a weak dependence on Sersic 
index to), the line-of-sight velocity dispersion varies as I? -1 / 2 , and the velocity anisotropy is (3 ~ —0.24 
to —0.28 depending on to. The excess stellar mass in the cusp scales approximately linearly with the 
black-hole mass, and is typically 0.5-0.85 times the black-hole mass. This process may strongly 
influence the structure of nuclear star clusters if they contain black holes. 


1. INTRODUCTION 


Most elliptical galaxies have surface-brightness profiles 
that are fit well by the empirical Sersic formula, given by 
equation (JT]) , with Sersic index to varying from 2 to much 
larger values (K ormendy et al.| 2009). In addition to the 
index m. the Sersic formula has just two parameters: In 


just two parameters: 
is the central surface brightness and R e is the effective 
radius, the projected radius containing half of the total 
luminosity. 

It has been observed ( Kormendy et al.||2009 ) that the 
Sersic profile fits elliptical galaxies over a wide range of 
radii, but there are often deviations from the fitting for¬ 
mula at small radii. These deviations are of two kinds: 
either “missing light” or “excess light” near the centers 
compared to the inward extrapolation of the Sersic pro¬ 
file. Missing-light or “core” ellipticals generally have 
higher luminosities {My 1$ —21.7) than excess-light el¬ 
lipticals {My > —21.5). Figures [l] and [2] taken from Ko¬ 
rmendy et al. (2009), show the surface-brightness profiles 


of an excess-light and a missing-light elliptical galaxy. 

The formation of a core is generally attributed to the 
dynamical excitation of stellar orbits due to a binary 
black hole (hereafter BH), although there is little or no 
direct evidence for this hypothesis. Excess light is gen¬ 
erally attributed to late star formation from gas that 
has become concentrated near the galaxy center. The 
division of ellipticals into these two classes, and the ex¬ 
planations given above for this dichotomy, date back at 
least to iFaber et al.l (119971) . 


Most elliptical galaxies appear to have massive BHs at 
their centers. In this work we assume that these BHs 
form slowly by the accretion of gas from larger radii. In 
this case the surrounding stellar system will contract as 
the central mass grows, because the stellar orbits con- 
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r'/ 4 (orcsec 1 / 4 ) 


Fig. 1 . — Surface-brightness profile of the dwarf galaxy NGC 
4459, an excess-light elliptical. The curve shows the best- 
fitting Sersic profile, with index m = 3.17 (|Kormendy et ah] 


2009). Reproduced by permission of the AAS. 


serve their adiabatic invariants or actions. If the galaxy 
began with a Sersic surface-brightness profile, this pro¬ 
cess will naturally produce excess light relative to the 
Sersic profile near the center of the galaxy, and the cor¬ 
responding excess stellar mass will be simply related to 
the mass of the BH and the properties of the initial Sersic 
profile. The main goal of the paper is to work out these 
relations and to compare the predictions of this simple 
model to the observations. 

Pioneering studies of this process were carried out by 


Peebles (1972) and Young (1980). These studies assumed 


that the galaxies were spherical, as do we, but in contrast 
they assumed that the stellar distribution function (DF) 
near the center of the galaxy was Maxwellian before the 
formation of the BH. They found that as the BH grew 
the stellar DF developed a cusp in which the density var¬ 
ied as p ~ r -3 / 2 , corresponding to a surface brightness 
I{R) ~ I? -1 / 2 . However, a Sersic profile does not have a 
Maxwellian DF, and the properties of the cusp formed by 
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the adiabatic growth of a central BH depend strongly on 
the DF (Quinlan et al(| |199 5). In general we expect (and 
shall find) the cusp formed in an initial Sersic model to 
be steeper than the one found by Peebles and Young. 



r '/ 4 (arcsec 1 / 4 ) 

Fig. 2.— Surface-brightness profile of the giant galaxy NGC 
4472, a missing-light elliptical. The curve shows the best¬ 
fitting Sersic profile, with index m = 5.99 (|Kormendy et ah] 


2009). Reproduced by permission of the AAS. 


In fj2] we review the properties of Sersic models and 
derive the dynamical quantities of interest (density, po¬ 
tential, DF, velocity moments) assuming spherical sym¬ 
metry, constant mass-to-light ratio Y, and isotropic ve¬ 
locity dispersion. We then slowly add a central point 
mass representing the BH and evolve the system un¬ 
der the assumption of adiabatic invariance of the stel¬ 
lar orbits. The numerical methods are described in (0 
The surface-brightness distribution, anisotropy parame¬ 
ter, excess mass, and line-of-sight velocity-dispersion dis¬ 
tribution of the adiabatic models are described in §(0and 
0 an d the observations in jj6j The conclusions are in 0 


The luminosity density for given surface brightness 
I(R) can be obtained by the Abel integral 


j(r) 


1 f°° d R d I 
tt Jr \/R 2 - r 2 d R 


(3) 


If we assume that the galaxy has constant mass-to-light 
ratio T, then its mass density p(r) = T j(r). Substituting 
equation 0 in equation 0 and introducing the dimen¬ 
sionless radius s = r/R e , we get 


Pm(s) = T^-pm(s) = ^Pmis) where E 0 = TJ 0 , 


Pm{s) = 


rl exp[— b(s/x) 1 ^ m ] dx 


irms( m 1 )/ m 


5 l/n 


t/l — 


(4) 


where x = s/rj. For m > 1 the dimensionless density 
Pm{s) diverges at the origin as 


Pro(s) ~ 2Tvrn a (™~ D/m B ib F 1 ~ V™)] (5) 

where B is the beta function. The dimensionless density 
Pm(s) is plotted in Figure 0 for several values of m. 



2. INITIAL SET UP 
2.1. Density 


The surface-brightness distribution of a spherical galaxy 
is modeled well by the Sersic law (|Sersic 1963 Sersic 
1968) |Ciottip99l| ) 


I(r)) = Iq exp(— bp 1 /™) where r] = R/R e (1) 


and m is the Sersic index. Here R e is the effective radius, 
the radius on the sky that contains half the total light. 
This definition of R e requires that the parameter b is 
given by 

2 7(2m, b) = T(2m) (2) 


where 7 is the incomplete gamma function. The function 
b(m ) is well fitted by the linear interpolation b(m) = 
2m — 0.324, with relative error less than 0.001 for 0.5 < 
to < 10 (Ciotti 1991). The exact values of b(m) for 
several values of to are given in Table [l] As to —* 00, 
b/m —> 2 and I{R) ~ R ~ 2 over a wide range of radii 
around R e . 


Fig. 3.— The dimensionless density p m (s) (eq.pj) for 2 < m < 11. 
Here s = r/R e , the ratio of the radius to the effective radius. 


2.2. Potential 

The potential is obtained by solving the Poisson equation 

V 2 $ = 4y- ( r2 7n) = AnG P- ( 6 ) 

r z dr \ dr j 


If we define c f > m (r) = 47rG T,oR e ^ m (s), then in dimen¬ 
sionless variables we have 


y d_ / 2 d 4» m A 
s 2 ds ds J 


— Pm • 


(7) 


In integral form, 


$m(«) 

M m (s ) 





„/2 


ds' 



where 


( 8 ) 
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is the mass of the system measured in units of 47 t Y,qR^. 
One cannot solve the integral in general for the density 
given b y equat ion d4|) but we can obtain the central po¬ 
tential (Ciotti 199ijb which is useful in what follows: 


<F(0) = —4GY / dRI(R) = — 4-7tG / r p(r) dr (9) 


from which the dimensionless central potential is found 
to be: 


$m(0) 


F(m + 1) 

nb m 


( 10 ) 


Thus the central potential is finite even though the den¬ 
sity diverges at the center for to > 1. The values of the 
dimensionless central potential and dimensionless mass 
for several values of the Sersic index m are given in Ta¬ 
ble on 

Poisson’s equation |7]) has to be solved numerically 
with the initial condition (lo| , taking d4> m /ds = 0 at 
s — 0. Plots of the potential for different values of to are 
shown in Figure [4] 



Fig. 4.— The normalized potential < t> m (s)/|<t> m (0)| (eqs. [s] and 
110| ) for 2 < m < 11. Here s = r/R e , the ratio of the radius to the 
effective radius. 


TABLE 1 _ 

The dimensionless parameters b, |<I> m (0)| and M m for 

DIFFERENT VALUES OF THE SERSIC INDEX m. 


m 

b 

l*m(0)| 


2 

3.67 

4.17x 10 -2 

2.91x 10 -2 

3 

5.67 

1.04x 10 -2 

5.41x 10 -3 

4 

7.66 

2.21x 10 -3 

8.42x 10 -4 

5 

9.66 

4.52x 10~ 4 

1.27x 10 -4 

7 

13.66 

1.80x 10 -5 

2.74x 10 -5 

9 

17.66 

6.88x 10 -7 

5.68x10-® 

11 

21.66 

2.57x 10 -8 

1.17X10 -9 


2.3. Distribution function 

We shall assume for simplicity that the initial phase- 
space DF of the stars in the galaxy is ergodic, that is, the 
velocity distribution is isotropic at every point in space 
and the DF depends only on the energy E. In such a 
case, the DF can be calculated from the known density 
and potential of the galaxy as follows: 


n } VSitzJe d $ 2 


( 11 ) 


Using the dimensionless potential and density, we can 
write the DF in dimensionless form as 


fm{E) 

fm(E) 


(47rG) 3 / 2 (E 0 -Re) 1/2 
1 f° d$ d 2 p m 

V8tt 2 Je - E d$ 2 


where 


( 12 ) 


The DF is plotted for various values of to in Figure [5] 
One can see that for all values of to the DF is positive, 
and that it diverges as E -A 4>(0); it is straightforward 
to show that 


fm(E) 


E-$ m ( 0) 


—a 


5 m + 1 

2 (to + 1)' 


(13) 



Fig. 5.— The DF / m (E) as a function of normalized energy 
E/&m( 0) for 2 < m < 11. 


2.4. Velocity moments 

For the remainder of this paper we will drop the tildes 
on E , <f>, p, etc., and it should be understood that these 
quantities are dimensionless. According to Jeans’s theo¬ 
rem, a stationary spherical DF can only depend on the 
energy E = \vf + + <I>(s) and the angular momentum 
L = svt, where v r and Vt are the radial and tangential 
velocities. We can obtain the density, and the radial and 
tangential velocity dispersions for this DF as 
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P(s) 

Vr(s) 

V?( s ) 


4n f° d E[2(E-$)]V 2 f d yyf^ ^y) t 
J$(s) Jo V 1_ 2/ 2 

f dE [2(E - $)] 3/2 f 1 dyyVl-y 2 f(E,L max y) : 
J<S>(s) JO 

dE [2 (E - $)] 3 / 2 f 1 dV y3 l^ ax V) , 

J<p(s) Jo y/1 - y 2 


An 

P(s) 

47T 

p(s) 


(14) 

(15) 

(16) 


where 

L mSLX (E,s) = {2s 2 [E-<S>(s)\} 1 / 2 , y = L/L maa t (17) 


Lmax is the maximum angular momentum for an orbit 
of energy E that passes through radius s. Note that 
the velocity dispersions in the 6 and </> directions are 
Vq = Vg = \v 2 . For numerical work, the inner integral 
over y can be evaluated using Gauss—Chebyshev quadra¬ 
ture and the outer integral over E by Gauss-Legendre 
quadrature. 

We then project the density and velocity moments on 
the plane of sky to obtain the surface density 


Efa) = 2 



s p{s) ds 
a/s 2 — if 


(18) 


and the line-of-sight velocity dispersion 


oj|(f?) 


2 r°° s p(s) ds 


The anisotropy parameter is 




(19) 


/3(s) = 1 - 



( 20 ) 


which is 0 for an isotropic velocity distribution, — oo for 
circular orbits and +1 for radial orbits. 



Fig. 6.— The initial isotropic velocity dispersion of the 
galaxy for varying Sersic index m. 


The initial Sersic models are ergodic, so the velocity 
distribution is isotropic, fi = 0, and the DF depends 
only on energy, not angular momentum. The velocity 
dispersion u 2 (s) = |u 2 (s) is shown in Figure [gJ 

2.5. Actions 

The radial and azimuthal actions are 
I r = 2 f + y]2[E- $(s)] - L 2 /s 2 ds I t = 2 n L, 

J S — 

( 21 ) 

where s_ and s+ are the pericenter and apocenter dis¬ 
tances of the orbit, given implicitly by 

2s 2 ± [E- $(s±)] - L 2 = 0. (22) 

Since actions are conserved during the adiabatic growth 
of the central BH, the DF is also conserved at each point 
in action space. In terms of energy and angular mo¬ 
mentum, a star that is initially at ( E , L) in phase space 
when the potential is < F(s) will move to a new phase 
point (E*,L) when the potential changes to $*(3). By 
equating the initial and final actions, we find the relation 
E = E(E*,L). Then the new DF f*(E*,L) is related to 
the original DF f(E, L) by f*(E*,L) = f[E(E*, L), L\. 

3. NUMERICAL METHODS 


The density, DF, potential and other galaxy properties 
are evaluated on a grid: pt = p(s,;), >!>; = <b(si), fjk = 
f(Ej,Lk), and so on. The radial grid is equally spaced in 
log(s) with roughly 200 points between s m i n = 10 -6 and 
Smax = 2 x 10 2 . The energy grid points Ej are chosen 
to match the potential T, at the radial grid points; since 
>I>(r) varies as the BH grows, the locations of the energy 
grid points Ej also vary. The grid points for angular mo¬ 
mentum are spaced uniformly between x = 0 and x = 1 
where x = L/L C (E) and L C (E) is the angular momentum 
of a circular orbit at the given energy. Typically we use 
200 grid points in energy and 20 in angular momentum. 

Following Young (1980), the program executes the fol¬ 
lowing steps for a given Sersic index m: 


1. The density pi = p(si) is obtained from equation 
0 at all points on the radial grid (recall that we 
have dropped the tildes on all dimensionless quan¬ 
tities) . 

2. The gravitational potential <fq = 4>(si) is obtained 
from Poisson’s equation (|8j) at all points on the 
grid. 
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3. The initial DF fjk = f(Ej ) is obtained from pi and 
<f>i using equation (12). The initial values fjk are 
independent of the angular momentum XkL c (Ej ) 
since the initial DF is ergodic. 


4. The radial action Ij k = I r (Ej,Xk) is calcu¬ 
lated for the potential $(s) at each point on the 
energy-angular momentum grid [Ej,Xk L c (Ej)], us¬ 
ing equation ( 21 ). 


5. A BH of mass M, is added at the center of the 
galaxy, thereby modifying the potential to <f>*(s) = 
<F(s) — M,/s. 


6. In this new potential we compute the angular mo¬ 
mentum of circular orbits as a function of energy, 
L* C {E*). 

7. A new energy-angular momentum grid 
[E*,Xk L* C (E*)} is defined, in which the new 
energies E* match the new potentials 

at the radial grid points. The new action 
I* k = I*(E*,Xk ) is computed on the new energy- 
angular momentum grid [E*, Xk L*(E*)]. 

8. Since actions are conserved, a star at (E, x) 

evolves during the growth of the BH to ( E*,x*) 
where I*(E*,x*) = I r (E,x), x* L*(E*) = 

xL c (E). We solve for the original (E,x) at each 
point on the new energy-angular momentum grid 
[E*,Xk L*(E*)} and then compute the new DF as 
f*(E*,x*L * c ) = f[E(E*,x*L* c )\. 

9. Using the DF f* and potential <f>* we compute the 
resulting density p** from equation ( |14| ) and the 
potential $** from Poisson’s equation (tST. We then 
replace p* and 4>* by p** and <&** and return to step 
6 . 


After about 20 iterations this algorithm converges, 
yielding a DF for the stellar system from which quan¬ 
tities of interest like density, surface brightness profile, 
and velocity dispersion can be computed. 

4. EXCESS MASS 

We wish to quantify the excess light or stellar mass accu¬ 
mulated near the center of the galaxy as the result of the 
adiabatic adjustment of the stellar orbits to the growth of 
the BH. Since the excess has appeared at the center only 
because of the rearrangement of the stars, the total lu¬ 
minosity of the galaxy is unchanged: the density goes up 
substantially close to the center, and decreases by a much 
smaller amount at large radii. The definition of “excess 
mass” is therefore somewhat arbitrary, but most reason¬ 
able definitions will give very similar values. We define 
the “cross-over radius” as the value of the dimensionless 
radius 77 = R/R e at which the initial and final surface- 
density profiles cross (see Figure [7]). Let S( 7 y) and £*( 77 ) 


be the dimensionless initial and final surface-density pro¬ 
files, then 77 c is defined by the equality T,*(ri c ) = Y,(rj c ); 




Fig. 7.— Typical potential and surface density of a galaxy 
with and without a BH. The Sersic index m = 4 and the 
black-hole mass relative to the galaxy mass is M,/Mg = 
1.877 x 10“ 3 . 


we have verified that this solution exists and is unique. 
The excess mass is defined as 

= \ [ ??[S*(? 7 ) - 11 ( 77 )] dr]. (23) 

The factor \ appears in the above equation so 
that the dimensional definition of the excess mass is 
M e = 2n J^ C [E*(R) — T,(R)]RdR (cf. eq. J 8 j)— in other 
words it arises because the volume element is 4-7ts 2 ds 
but the area element is 27777 ^ 77 . 

5. RESULTS AND DISCUSSION 
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R/R e r/R e 


Fig. 8.— Surface-brightness profile (left) and anisotropy parameter (right) after the adiabatic growth of a BH in a galaxy with 
Sersic index m = 2.3, m = 3, m = 6, and m = 9. 
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The surface-brightness profile I(R) and anisotropy pa¬ 
rameter /3(r) are plotted in Figure^] for Sersic index m = 
2.3, 3, 6 and 9 for dimensionless BH masses varying from 
10 -5 to 5 x 10 -3 . The BH masses are chosen such that 
they are 10~ 4 to few times 10~ 2 of the galaxy mass (-Mg). 

5.1. Surface-density cusp 


The prominence of the surface-brightness or surface- 
density cusps shown in Figure [ 8 ] varies with the Sersic 
index to. For smaller m the cusp is more prominent, 
partly because the initial surface-brightness profile near 
the center is flatter. For large to, even with the largest 
BH masses, the cusp caused by the BH is hardly notice¬ 
able (see to = 9 case in Figure [ 8 ]). 

One can derive the approximate behavior of the 
surface-density cusp analytically ( |Quinlan et al. 1995 
Gondolo &; Silk|[l999 ) 


-r 


R 


l-r 


r = 2 


4 — y 


7 to + 2 
3 to + 1 


(24) 

where 7 = (to — 1 )/m is the exponent of the initial cusp 
(see eq. [ 5 J. For the Sersic index to = 2.3, p* ~ r ~ 2 29 
and £* ~ R~ 129 and for to = 9, p* ~ r " 2 ' 32 and £* ~ 
i? -1 - 32 . We see that the exponent of the central surface- 
density cusp does not vary strongly with Sersic index. 
Thus the cusp becomes less prominent as m increases 
and the initial profile becomes more cuspy. 


5.2. Anisotropy and velocity cusp 


The line-of-sight velocity dispersion (eq. 19) is shown in 
Figure [ 9 J Both the radial and tangential dispersions v 2 
and v 2 scale as r _1 near the BH and this scaling is evident 
in Figure [9] In contrast, the velocity dispersion goes to 
zero at the center for all values of the Sersic index to for 
a galaxy without a central BH. 

The anisotropy parameter f} is plotted for to = 2.3 
to to = 9 in Figure [8j The plots show that as a re¬ 
sult of the BH growth the stellar DF develops a tan¬ 
gential anisotropy — /3 ~ 0.24-0.28 near the center. The 
anisotropy depends only weakly on the Sersic index, and 
is independent of the BH mass: as the BH mass grows 
the anisotropy simply extends to larger distances. 

Using an approximate expression for the radial action 
near t he center of a Sersic model due to lGondolo fe Silk 
(1999), we compute in the Appendix the ratio of the 
radial and tangential velocity dispersions and thus the 
anisotropy parameter f) (eq. 201. When sufficiently close 
to the center, the anisotropy ratio is independent of the 
BH mass and the radius, and is determined solely by 
the initial Sersic index to. The values of /3 for m = 2.3 
and to = 9 obtained from equation (A6) are —0.313 and 
—0.306, respectively, and the corresponding values from 
Figure[8]are —0.285 and —0.247. The constant tangential 
anisotropy near the center contrasts with the situation 
for initial galaxies having analytical cores as described 


Goodman & Binney (1984) and Quinlan et al. (1995): 

0 


with analytic cores the initial DF is constant as r 
so the ratio v 2 /v 2 tends to \ and fd tends to zero. 



R/R, 

Fig. 9.-- Line-of-sight velocity dispersion after the adiabatic 
growth of a BH in a galaxy with Sersic index m = 2.3, m = 4, 
m = 6, m = 9. The dashed line is for the initial model 
without a BH. 
























































mass shown in Figure [TT] can be fitted to a cubic polyno¬ 
mial, 


V = 


where 


Me 

M. 


ax 3 + t>x 2 +cx + d ; 


X = log 


M. 

Mg 


(25) 


where the coefficients a, b, c and d are functions of 
m alone. Table [2] gives the fitted values of the coeffi¬ 
cients of the cubic polynomial; the ratio of black hole 
mass to the galaxy mass at the maximum of M e /M. 
(M™ ax /Mcj); and the maximum deviation of the fitting 
function ya t from the numerically calculated values j/num, 
that is, max(| 2 / fit - y num \). 


6. COMPARISON TO OBSERVATIONS 


Fig. 10.— Excess mass vs. BH mass, both scaled to the 
mass of the galaxy. The dashed line is M, = M e . 



Fig. 11. — Ratio of the excess mass to the BH mass, as a 
function of BH mass. The abscissa is plotted on a log scale. 


We now ask whether our results are consistent with ob- 



Fig. 12.— Surface-brightness profile of the galaxy with Sersic 
index m = 2.3, plotted with respect to ( R/R e ) 1/1 . The dashed 
line is M, = 0. 


5.3. Excess mass 

In Figure 10 we plot the excess mass (as defined in 
with respect to the BH mass. We observe that the excess 
mass scales approximately linearly with the BH mass, 
although with some offset and curvature in the relation. 
To make these features more apparent we plot in Figure 
EH the ratio of excess mass to BH mass as a function of 
the BH mass. 

We see that the excess mass M e is always smaller than 
the BH mass M, for all values of the BH mass and all 
values of the Sersic index m. There is a peak in the 
excess mass as a function of BH mass for all values of 
the Sersic index, which occurs when the mass of the BH 
is between about 10~ 3 and 10~ 15 times the mass of the 
galaxy. The maximum of the ratio M e /M, increases as m 
decreases, consistent with the stronger visual appearance 
of the cusp in Figure [8] at smaller Sersic index, but this 
maximum varies by less than a factor of two over a wide 
range of Sersic indices. 

The dependence of the ratio of the excess mass to BH 



at the lower end of the observed range. This result is 
consistent with the top panel of Figure 42 of |Kormendy| 


et al.l (2009), which shows that M e /M . ranges between 


1 and 50, with typical values around 10. In other words 
the adiabatic model correctly predicts that the excess 
mass should be roughly proportional to the galaxy or 
BH mass, but the predicted magnitude is too small by a 
factor 5-10. We conclude that the excess light in ellip¬ 
ticals does not arise from adiabatic growth of a central 
BH. This conclusion is consistent with the observation 
( |Kormendy et ah 2009) that the excess light often has 
disky isophotes, suggesting that it formed via dissipative 
processes rather than by adiabatic ones—most likely late 
star formation due to gas that has collected in the central 
region of the galaxy (e.g., Hopkins et al. 2009). 
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Nuclei or nuclear star clusters are a different form of 
excess light at the centers of galaxies. Kormendy et 
al. (Section 9.7) make a distinction between the “excess 
light” in elliptical galaxies and the “nuclei” in spheroidal 
galaxies. In Figure [12] the surface brightness profile for 
a galaxy with Sersic index m = 2.3 is plotted using the 
same ( R/R e ) 1 ^ 4 abscissa as in Kormendy et al. Compar¬ 
ing this plot to Figures 25, 26, and 28 of that paper, we 
notice that the slope of the surface-brightness profile as¬ 
sociated with nuclei in spheroidal galaxies, which is much 
steeper than the profile associated with excess light in 
ellipticals, is similar to what we obtain in the adiabatic- 
growth model. Moreover, the median light fraction of 
nuclei in spheroidal galaxies is 0.3% (Kormendy et al. 


2009), much smaller than the fractional excess light in 


The observational relation between BHs and nuclear 
star clusters remains unclear. Cote et al. (2006) and 


Wehner & Harris (2006) have argued that nuclear star 


clusters and central BHs define a single smooth relation¬ 
ship between mass (of the cluster or BH) and galaxy 
luminosity, with clusters dominating at low luminosities 
and BHs at high luminosities. Kormendy et al. (20091 


argue that this correlation is accidental, but our under¬ 
standing of the correlation between BH masses and their 
host galaxy properties, particularly in low-luminosity 
and late-type galaxies, is too limited for a definite con¬ 
clusion. 

One attractive alternative to the formation of nuclear 
star clusters by the adiabatic growth of a BH is that 
the clusters form through in-spiral of globular clusters 


ellipticals and consistent with our model if the nuclei con- (Tremaine et al. 1975 Antonini 2013 Gnedin et al. 2014). 


tain typical BHs with mass M,/Mg — 0.5% (Kormendy 
¥Hofl20T3 |. 


In this case the BH could form before or after the bulk 
of the cluster formation. 


TABLE 2 

The fitted values of the coefficients a, b, C d IN equation (|25|) AND M” ax /M G , THE BH MASS at WHICH M e /M. IS LARGEST. 


m 

—a 

-b 

— c 

d 

M” ax /M G 

max (| 7/fit 2/num|) 

2 

1.21x 10 -2 

1.84x 10 _1 

0.77 

-0.140 

5.21 X 10~ 2 

6.34 x 10“ 3 

3 

5.25x 10 -3 

9.33x 10 -2 

0.46 

0.0598 

3.03 x 10~ 2 

5.183 x 10 -3 

4 

3.23x 10 -3 

6.36X10" 2 

0.356 

0.0897 

1.76 x 10~ 2 

3.107 x 10 -3 

5 

2.47x 10 -3 

5.13x 10 -2 

0.312 

0.0784 

1.09 x 10~ 2 

1.115 x 10~ 3 

6 

1.88x 10 -3 

4.14x 10 -2 

0.272 

0.0810 

6.78 x 10~ 3 

0.568 x 10" 3 

7 

1.39x 10 -3 

3.31xl0 -2 

0.236 

0.0934 

4.21 x 10~ 3 

0.761 x 10“ 3 

8 

1.03x 10 -3 

2.64x 10 -2 

0.205 

0.1070 

2.65 x 10~ 3 

0.968 x 10" 3 

9 

9.19x 10 -4 

2.41xl0 -2 

0.194 

0.0970 

1.83 x 10~ 3 

0.574 x 10~ 3 


7. CONCLUSIONS 


We have studied the effect of the adiabatic growth of a 
central black hole on the R}/ m or Sersic model of spher- 
ic al ga laxies, following the methods described by |Young| 
(1980). The black hole induces a surface-brightness cusp 


at the center of the galaxy, which can be described as 
“excess light” or “excess mass” above the inward extrap¬ 
olation of the best-fit Sersic profile for the outer galaxy. 
At the smallest radii the surface brightness is found to 
vary with radius as /(A) ~ A -1 ’ 28 to /(A) ~ A” 132 
for to varying from 2 to 9. Increasing the BH mass in¬ 
creases the strength of the cusp without changing the log¬ 
arithmic slope near the center. The line-of-sight velocity 
dispersion is found to scale as A -1 / 2 , which is the ex¬ 
pected behavior in a Keplerian potential; the anisotropy 
parameter /? near the center is between —0.24 and —0.28 
(tangential anisotropy). 

We calculated the excess mass (defined as the mass 
interior to the radius on the sky where the initial and 


final surface densities were the same) for Sersic models 
with varying indices to. The excess mass is generally 
between 0.4 and 0.85 times the black-hole mass. 

If the typical black-hole mass is ~ 0.5% of the stellar 
mass in the galaxy, the excess mass produced by adia¬ 
batic contraction is ~ 5-10 times smaller than the excess 
mass determined from Sersic fits to the photometry of el¬ 
liptical galaxies, but similar to the masses of nuclear star 
clusters in low-luminosity galaxies of all Hubble types. If 
black holes form by slow accretion of gas then adiabatic 
contraction may play an important role in determining 
the properties of nuclear star clusters. 
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APPENDIX 


The goal of this Appendix is to determine the velocity anisotropy near the center of the cusp formed by the adiabatic 
growth of a central BH. 

The density for a spherically symmetric Sersic profile of index to can be approximated near the center as (eq. |5j) : 

Pm{s) = p 0 s -T where 7 = (m — 1 )/to, 

and the corresponding potential is 


$m(s) - $ m (0) = 


Po 


S 2 7 = 4>o S 


2-7 


(3- 7 )(2-7) 

The phase-space DF for the above density and potential pair is given by equation (13): 

/„(£)-[E-$ m( 0)r», « = = 

An approximate expression for the radial action (Gondolo & Silk 1999) accurate to 8 % for all to > 1 is 


(Al) 

(A2) 


Ir{E,L) — —T \ - -T + \/2 <^) 0 


£-$m(0) 


(4- 7 )/(4-2 7 ) 


(A3) 


(A4) 


where 


A = 


2 i /(2- 7 )( 2 _ 7)1/2 

(4 _ 'y)(4— 7 )/(4—2 7 ) : 


d = 


. (2-7) 

B[ 1/(2- 7 ), §]' 


Once the black hole grows adiabatically, the galaxy potential at the center is approximately Keplerian. The radial action 
is now given by I*(E*,L) = 2tt(M,/\/— 2 E* — L). The final DF f*(E*, L) is obtained by solving I r (E, L) = I*(E*,L ) 
for E as a function of E* and L. Using the equation E* = \v 2 — M t /s and L = svt to eliminate E* and L in favor 
of the total and tangential speeds v and v t , we find that the dependence of the final DF on velocity at a given radius 
can be written as 

1 \ T 2 

— — 1 \x sin %p + 


where 


S = 


r 

6—7 
4—7 


A d 


1 -5 


2V x ‘L-x 2 _ 

= \/ 2M, s ; vt = v sin ijj . 


(A5) 


The quantities x and x m are defined following Goodman & Binney (1984). 

Multiplying f* by v 2 and v 2 and integrating over all velocities we obtain the ratio of velocity dispersions and the 
anisotropy parameter /3: 

n -8 


d y y 


dip sin ip cos 2 ip 


1 

A d 


— 1 I y sin ip + 


1 


2vT^ 


y 


■J 2 2(1 - P) 


d y y 


d 'ip sin 3 ^ 


— - 1 ) y sin ip + 




y- 


-8 


(A 6 ) 


where y = x/x m . The only approximation in deriving this result is the approximate form for the radial action, equation 
(A4). The quadratures are straightforward to evaluate numerically. 




































